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Abstract 

In this work, we evaluate the lifetimes of the doubly charmed baryons Sj,, and f^. We care- 
fully calculate the non-spectator contributions at the quark level where the Cabibbo-suppressed di- 
agrams are also included. The hadronic matrix elements are evaluated in the simple non-relativistic 
harmonic oscillator model. Our numerical results are generally consistent with that obtained by 
other authors who used the diquark model. However, all the theoretical predictions on the lifetimes 
are one order larger than the upper limit set by the recent SELEX measurement. This discrepancy 
would be clarified by the future experiment, if more accurate experiment still confirms the value 
of the SELEX collaboration, there must be some unknown mechanism to be explored. 
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I. INTRODUCTION 



The quite large difference of the lifetimes between and D° and the lifetimes close 
to each other for B^ and B° are well explained by taking into account the non-spectator 
effects l|. This success implies that the mechanism which governs the reactions at quark 
level is well understood. When we apply the mechanism to the heavy baryon case, some 
problems emerge. The famous puzzle in the heavy-flavor field that the lifetime of is 
remarkably shorter than that of B meson is much alleviated recently when the operators 
of higher dimensions are taken into account j^, 3|. The more recent experimental value of 
the ratio t(Ai,)/t(B ) = 1.041 ± 0.057 4] is close to the theoretical evaluation^. However, 
in the theoretical works, one can notice that the evaluation of hadronic matrix elements is 
still very rough and based on some approximations. The possible errors brought up by the 
uncertainties in the hadronic matrix elements are still uncontrollable. In our recent work {J, 
we find that the short-distance contributions to the branching ratio of A5 — > A7 which is 
evaluated in the PQCD approach, are much smaller than that from long-distance effects. 
Therefore, even though one has a full reason to believe that the low-energy QCD should 
solve the discrepancy if it exists, he must find a proper way to deal with the hadronic matrix 
elements. 

The observation of doubly charmed baryon by the SELEX Collaboration at 
FERMIL AB ] provides an opportunity to investigate the hidden problems. Hopefully the 
study may shed some lights on the unknown non-perturbative QCD effects which result in 
obvious difference between baryons and mesons. Because contains two heavy quarks, 
by the heavy quark effective theory (HQET) the situation may become relatively simple 
and clear compared to the case of A& or A c which possesses only one heavy quark. Thus a 
careful study on the is necessary and interesting. Several groups already investigated 
the two-heavy-flavor baryons a long time ago[3, Isj]. In their work, the evaluation of the 
hadronic matrix elements is based on the quark-diquark structure of the baryons. This is 
definitely reasonable, it is believed that two heavy quarks can constitute a more stable and 

. However, since charm quark, even b-quark, is not so 
heavy that the degree of freedom of the light flavor can be ignored, the diquark scenario 
may bring up certain errors, especially when evaluating lifetimes of baryons, because only 
inclusive processes are concerned. In this work, we do not use the diquark picture, but 
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instead, adopt a simpler non-relativistic model for the baryon and re-evaluate the hadronic 
matrix elements. As a by-product, one can compare the results by the diquark picture with 
that by the three valence-quark picture. It may help us to better understand the diquark 
picture and its application range. 

The advantage is obvious, that we only concern the inclusive processes in terms of the 
optical theorem when calculating the lifetime. Therefore, we do not need to deal with 
the hadronization to light hadrons. The only non-perturbative effects come from the wave 
function of the heavy baryon. Moreover, since there are two heavy quarks in the baryon, 
the relativistic effects are not so significant and the framework of non-relativistic harmonic 
oscillator model might lead to a reasonable result. 

Moreover, at the quark level, we carry out similar calculations as that in the literature, 
but we keep some new operators which are CKM suppressed and contribute to the lifetime. 
They appear at the non-spectator scattering at order of in heavy quark expansion(HQE). 
Later, our numerical results show that their contributions are indeed very tiny to make any 
substantial contributions. 

All the concerned parameters in the model are obtained by fitting data, therefore we 
avoid some theoretical uncertainties and obtain reasonable results. Comparing these results 
with data, we may gain information about the the whole picture. 

Our paper is organized as follows. In Section. II we derive the formulation for the lifetimes 
of S+, 5+ + and fl^ c which include the non-spectator effects. In Section. Ill, we use a simple 
model, i.e. the harmonic oscillator, to estimate the hadronic matrix elements. In Section. IV 
we present our numerical results along with the values of all the input parameters. The last 
section is devoted to our conclusion and discussion. 

II. FORMULATION FOR LIFETIMES OF H+, H++, tt+ 
A. Spectator Contribution to Lifetimes of S+, fi+ 

The lifetime is determined by the inclusive decays. Thus one can use the optical theorem 
to obtain the total width (lifetime) of the heavy hadron by calculating the absorptive part 
of the forward-scattering amplitude. 
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The total width is then written as 



T(H Q -> X) 



rriHr 



Im / d 4 x(H Q \f\H Q ) 



where 



2m 



■(H Q \t\H Q ), 



(1) 



f = T{i£ eff (x),C eff (0)} 



(2) 



and £ e ff is the relevant effective Lagrangian. I/tuq is the expansion parameter, and the 
non-local operator T is expanded as a sum of local operators and the corresponding Wilson 
coefficients include terms with increasing powers of 1/rriQ. Definitely, the lowest dimen- 
sional term dominates in the limit nig — > oo and it is the dimension-three operator cc. 
The total width of a charmed hadron H c is determined by lm(H c \T\H c ) 
normalization 



X 01 ] with a proper 



r(H c 



f) 



G\m\ 

1927T 3 

+<*(/) 



VfciCM| 2 {c 3 (/)(^c|cc|^c) 

(H c \Ba^G^c\H c ) 



mi 



£<£></) 



(H c |(cr jg )(gr (C )|/f c ) 



O(-i)}, 



(3) 



where the coefficients depend on the masses of the internal quarks in the loop. The 
coefficient c 3 (/) has been calculated to one- loop order Q, Q, Q whereas the coefficient 



05(f) is evaluated at the tree level 15|, [l6|. Vckm is the Cabibbo-Kabayashi-Maskawa mixing 
matrix elements and G M „ is the gluonic field strength tensor. Since the third term involves 
light quarks, it can be different for charmed hadrons with various light flavors. Thus, the 
difference appears at the \jm\ order and in the hadronic matrix elements of four-quark 
operators. The contributions at orders higher than \jm\ are neglected. 

To the lowest order, the main contribution comes from the heavy quark(charm quark) 
decays, while the light flavors are treated as spectators. The contributions are due to the 
semileptonic and the nonleptonic decays as follows: 



Tic 



c— >slv ^ ^ ^ c—*sqq' 

=e,/x q(q')=u,d,s 



(4) 



The semileptonic and nonleptonic decay rates of the c quark up to order l/m^ has been 
evaluated by many authors[171|. and here we would directly use their results. 
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B. Non-spectator Contributions to Inclusive Decays of S+, 0+ 

The total width of hadrons which involve at least one charm quark c can be decomposed 
into two parts 

J^^^Jq > y) -^spectator _|_ -^nonspectator 

For the spectator scenario, the contribution to the total width of the (cccQ-baryon ground 
state S+, the (ccu)-baryon ground state 2+ + and the (ccs)-baryon ground state Vt+ C should 
be a sum of decays rates of two c— quarks individuallynamely 

r^ ec ~2rf ec , q = u,d,s. (6) 

To derive the non-spectator contributions for decays of E+, S+ + and we need the 



relevant effective Lagrangian: 



18| 



AG 

^e// (A* = m c) = — ^{^^[C r i(Ai)ay i Lc«7 M Ld+C2(/x)M7 M Lcs7 /1 Ld] 
+V cd V* d [C 1 (fi)d^Lcu lll Ld + C 2 (n)u^Lcd^Ld] 
+V C sV* s [C 1 {fi)s^Lcuj tl Ls + C 2 {n)u^Lcs^ IJL Ls] 
+V CS sry v M.'fLl} + h.c. (7) 

l=e,fi 

where L denotes — ^p-. 

(i) The inclusive decays of S+: 

There are four diagrams which contribute to the the width of S+, as shown in Fig.l. Here 
we also include the Feynman diagrams which are CKM suppressed. Fig l.(a),(c) are the 
W-exchange diagrams (WE), while Fig l.(b),(d) are the pauli- interference diagrams (PI). 
Here Fig l.(d) is arisen from the semi-leptonic decay of the charm quarks with the d— quark 
in S+. For the WE-type diagrams, we derive the contribution to the width as 

„ „+ 2G 2 

^WE = ~ (1 Vcs 1 | Kid I C(z s+ , Z u+ ) + | Vcd\ \V u< i\ C(z u+ , Zd+))P + 

{[C 2 M + Oi^crfLcd^Ld + 2C 1 ( f i)C 2 (fi)cYLdd^Lc}, (8) 

2 

where P + = p c + pd, z q+ = pf (<? = u,d, s). The definition of the function C(zi, z 2 ) is 

C{z u zz) = ~[-2(x 3 2 - xl) - {x 2 2 - x?)(3 + 2z x - 2z 2 ) + 4^(x 2 - Xl )], (9) 
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where x\ t2 = z 2)t\ / (^+zi-z 2 ) 4^ ^ ^ e expressions q and q are free field opearotors of 

quark and antiquark, and we will show in next section that all the non-perturbative QCD 
effects are included in the wavefunctions. Their explicit expressions are given as 

d 3 k 



q = J WfE q ^ K(kK(k)e-^ + d+ a (kX(k)e***) (10) 
* = I T^TT E (4(W) eita + <(kK(k)e- ikx ) ■ ( n ) 

For S+ , g=c, it. 

The contributions from the Pauli-interference(PI) non-spectator diagrams to the width 
of S+ are: 

+ 2f 2 

fp? = -^{|Kd| 2 |^| 2 P^«-,^-)[^ 

+2C 1 (/x)C 2 (/x)c7^L ?J 7 l 'Lc] + 2|K d | 2 F^(0, zx-)^ Ldd-f Lc} , (12) 

2 

where = pr(q = u,d, e, /i) and P_ = p c — Pd- The definition of the function F^zi, z 2 ) is 

3 

P^i, z 2 ) = -[2(x 3 2 - x?) - -(2 + z! - ^ 2 )(x 2 - x 2 ) + 3(x 2 - xi)]P 2 ^ 

+ [2 ( X 3 _ x 3 } _ 3{x 2 _ ^ (13) 

where the definitions of z\ and z 2 are the same as before. 

(ii) The inclusive decays of 

The non-spectator contribution to the width of H+ + come from the diagrams shown in Fig.2. 
That is caused by an interference of the produced u— quark from decay of one of the charm 
quarks with the u— quark in Here we also include the CKM suppressed Feynman 

diagrams. The contribution is 

++ 2Gr 2 

Tpj? = ~ {l^csl I 'Kid I Fuu[Zs—i z d—} ~\~ \^cs\ |Kis| Fiu/\Zs— > ^s— ) 
D7T 

+ |Kd| 2 |Kid| 2 

{z d _,z d _)} 

{C\{y)^Lv?u i 1 v Lc i + NCK^u^LccYLu + 2C x ^)C 2 {[i)u^LccUu}, 

(14) 

where z_ = ^-(q = s,d), P_ = p c - p u . 

(iii) For the inclusive decays of 

The non-spectator contributions for f2+ not only come from the Pauli interference of the 
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s— quark produced in the non-leptonic, but also from the semi-leptonic decay of the charm 



quarks with the s— quark in the later one is suggested by Voloshin et al.[l9[. As above, 
here we include the CKM suppressed WE non-spectator diagrams. The WE non-spectator 
contribution to the width Q~t is 



A + 2G 2 9 

^WE = |Kis| \Vcs\ C(z u+ ,Z S+ )P_ 



{[Cf(fi) + Cl(fx)}cYLcs^Ls + 2CMC 2 {^)crLs- Slil Lc}, 

(15) 



where z q+ = -p4, q = u,s and P + = p c + p s . 



+ 



The PI non-spectator contribution to the width of f2+ is 

- n + 2G 2 F 2 2 2 2 

Tpr = 7, { I ^cs | l^udl Fny(yZ u —^ Zd—^j ~\~ \ Vcs\ \Vus\ Favi^u—T^s—)) 

{NC^^c^LssYLc + CKny^Ls^i^Ld + 2C 1 {fi)C 2 (fx)crf' x LssYLc} 
2C 2 

-2— -iI|y cs | 2 F^(0, Zl _)cYLssYLc, (16) 

where z 9 _ = -p4, q = u,d,s,e,/i and P_ = p c — p s . Sandwiching the operators between 
initial and final S+, Q+ states, we obtain the hadronic matrix elements: 

^WE/PI = (^cefi* = 0' S ) \ ^ WE /Pll^tcfi* = 0' s )) 

rpf = <S£ + (P = 0, S )|fg + |H++(P = 0,,)) 

T we/pi = (^e( p = 0,a)|f j^ /PJ |n+(P = 0,8)). (17) 



III. THE HADRONIC MATRIX ELEMENTS 

Because the hadronic matrix elements are fully determined by the non-perturbative 
QCD effects which cannot be reliably evaluated at present yet, we need to invoke concrete 
phenomenological models to carry out the computations. In this work, we adopt a simple 
non-relativistic model, i.e. the harmonic oscillator [2^]. This model has been widely em- 



ployed in similar researches 2l|, |22|, |23|, [2J, |25|, |26j. In fact, an advantage of the calculations 



of the lifetimes of heavy hadrons is that one does not need to deal with the hadronization 
process of lighter products (quarks or even gluons) and the heavy hadrons can be well 
described by such simple non-relativistic models, and the results are relatively reliable than 
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for light hadron decays. 



(i)The inclusive decays of H+: 
In the harmonic oscillator model, the wavefunction of 5+ is expressed as |S+) and 

I — 'cc(^*> ^0) ^ ~X_ spin, flavor^ color 

color, spin 

d 3 p P d 3 p\ty E +{p p , p A )|c i (p 9l , 

S qi ) 5 C j (P?2 ' S q2 ) 5 die (Pq 3 , Sq 3 ) ) • 

(18) 

The normalization condition for |5+(P, s)) is 

(£+(P, S )|H+(P', S ')> = (2tt) 3 ^5 3 (P-P')V, (19) 

where Xspin, flavor^ 'color are the spin-flavor and color wavefunctions respectively Their explicit 
expressions are 

X s =\,fiavor = ^( 2 l c T c T rf l) - l c t c l rf T> - l c l c t rf t>) (2°) 

color = ^ijk- (21) 

As is the normalization constant. The spatial wavefunction \E f H + is a three-body harmonic 
oscillator wavefunction and expressed as 

2 2 

= exp(-^- - ^\). (22) 
PV 2a 2 p 2a\ ! K ' 

Here a p and a\ parameters reflecting the non-perturbative effects. In the above expressions, 
the Jacobi transformations of pi, P2, P3 which are the momenta of the three valence quarks 
ccd, and variables p p , Pa, P are 

Pi — P2 Pi + P2 — P3 

P P = 7=— ,Pa = , = — ,P = Pl + P2 + P3- (23) 

We choose the center-of-mass frame of S+, i.e. (P=0) to calculate the hadronic matrix 
elements. Substituting the four-quark operators into the expressions, we obtain the non- 
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spectator WE contributions to the width of as 



r*~ = 6Air 2 G 2 F Pl(\V cs \ 2 \V ud \ 2 C(z s+ ,z u+ ) + \V cd \ 2 \V ud \ 2 C(z u+ , z d+ ))(CM - C 2 (/x)) : 



spin 



(24) 



and the PI contribution is 

r p/ = -^Gl^Vj^VjfF^Zu-, z d J)[-NC 2 {ii) + C 2 (/i) - 2C 1 (n)C 2 (n)} 
-2\V cd \ 2 F, u (0, Zl _)}\A B \ 2 [2(l + ^)f 2 J2 [ d 3 p p d 3 p x d% 



spin 



eXph 2^ - 2^ ]eXph 2^ 2a] ]Uc7 ^ d ^ 7 LUc ' 

(25) 

where the sum over spin means a sum over the polarizations of the three valence quarks of 
S+ with their corresponding C-G coefficients in the spin-flavor wavefunction. u q , u q denote 
the Dirac spinors of free quarks q and the expression is 



2m q I Q--P 

E q +m q 



2m q X \ J { ' 

in our case q denotes c and d quarks. 
(ii)The inclusive decays of 

The contribution from the PI non-spectator diagrams to the width of H+ + is 

r S c+ = -^7r 2 ^{|^ cs | 2 |Kd| 2 ^(^-,^-) + |V; S | 2 |K S | 2 ^(^-,^-) 

+ \V cd \ 2 \V ud \ 2 F flu (z d -,z d „)}(C 2 1 (n) - NC 2 M - 2C l ( f i)CM)\A B \ 2 [2(l + ^)] 3 / 2 



E/^W 3 P;ex P [-^-^]exp[-| 



2a2 2a? J 

spin " P * P A 

u c YLu u u u j u Lu c . (28) 
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Similar to the case of , the sum over spin means a sum of the polarizations of the three 
valence quarks of H+ + with their C-G coefficients. One only needs to replace u by d in p p , 
Pa and other expressions are similar to that for S+. 
(iii)The inclusive decays of 

The contribution from the W-boson exchange (WE) non-spectator diagrams to the width of 

nt is 



m s J 



spin 



exp[-— - ^|]ex P [-^ ^ ]« cTm L^ s7 ^ s , 

(29) 

whereas that from the Pauli-interference(PI) non-spectator diagrams is 

[-TVC^/i) + C 2 2 (/i) - 2C 1 (//)C 2 (/x)] - 2|V CS | 2 F M1/ (0, ^_)}|A B | 2 [2(1 + ^)f 2 

Tfle 



£ / rf 3 Pp rf 3 p A rf 3 P ;exp[-^ - ^]exp 



P* pL r Pi + V 1 + ^ - P^) 



2(3^ 2a? 

sjnn P a P A 

u c j**Lu s u s YLu c . (30) 



The sum over polarizations is similar to that for 5+ and 



IV. INPUT PARAMETERS AND NUMERICAL RESULTS 

To obtain the decay amplitudes, we adopt the input parameters as follows [tI. I^: Gp = 
1.166xlO- 5 GeV~ 2 , \V CS \ = 0.9737, \V ud \ = 0.9745, Ci(m c ) = 1.3, C 2 (m c ) = -0.57, m c = 1.60 
GeV, m s = 0.45 GeV, m u = m d = 0.3 GeV, m* = 0.2GeV, ml = m* = 0, M =+ = M =++ = 
3.519 GeV, M 0+ = 3.578 GeV, M =+ * - M =+ = M =++ * - M =++ = M 0+ , - M 0+ = 0.132 

Si cc 1 — 'CC 1 — 'cc '—'cc ' — 'CC "CC "CC 

GeV. Here m q * denotes the current quark mass of flavor q. 

The non-perturbative parameters a p , a\ in the harmonic oscillator wavefunctions are 
selected as follows: for J/i}j, in ref. 20J], a 2 = 0.33GeV 2 , for mesons, a 2 = 0.25GeV 2 . For 
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TABLE I: The numerical results about the contributions from the different components and the 
evaluated lifetime for the doubly charmed baryons. For a comparison, in the following table, we 
list the corresponding lifetimes predicted by the authors of ref. [3] where the diquark picture was 
employed. It is noted that in ref.pj], the authors used various input parameters and obtained 
slightly diverse results, we take average values of the numbers in the table. There is only one 
datum for the lifetimes on r~+ given by the SELEX collaboration which is also listed the table. 



I— 'cc 


r spec (10- 12 GeV) 


rSio- 13 GeV) 


lti(10- 15 GeV) 


r s+(P s ) 


T~+ (ps) in ref. [7j 


exp(ps) 




2.01 


6.43 


-3.36 


0.25 


0.19 


0.033 


C++ 


r spec (10- 12 GeV) 




lt/ n (10- 12 GeV) 


(PS) 


(ps) in ref. [71 






2.01 




-1.02 


0.67 


0.52 






r spec (10- 12 GeV) 


r^(10" 14 GeV) 


lti(10- 12 GeV) 


r n+(P s ) 


T n+(P S ) ^ ref. [7] 






2.01 


4.25 


1.10 


0.21 


0.22 





the doubly charmed baryons, because a p reflects the coupling between two charm quarks, 
we set it to be the same as that for J/ip. a\ reflects the coupling of the light quark with 
these two charm quark, thus we can reasonably set it to be the same as a p in D-mesons. 

With these parameters as input, the lifetimes of the doubly charmed baryons can be 
evaluated out (see TABLE. I), if the non-spectator effects are taken into account. 



V. CONCLUSION AND DISCUSSION 

In this work, we evaluate the lifetimes of doubly charmed baryons with the non-spectator 
effects being properly taken into account. As argued in the introduction, to evaluate the 
lifetimes (the total widths), only the inclusive processes are concerned, and then the non- 
perturbative effects are all from the wavefunctions of the doubly charmed baryons. Due 
to existence of the two heavy charm quarks, the non-relativistic harmonic oscillator model 
should apply in this case. Mainly, we carefully calculate the contribution of non-perturbative 
effects to the lifetimes in the model, which are closely related to the bound states of the 
baryons. 

Our numerical results indicate that the non-spectator contributions to the lifetimes of 
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, H+ + and fi+ are substantial. The non-spectator contributions to the width of are 
mainly from the WE diagrams (the PI diagrams which contribute are CKM suppressed), 
since the WE contribution is constructive, therefore the lifetime of is much suppressed. 
By contraries, for E+ + and f2+, the non-spectator contributions are mainly from the PI 
diagrams and the net effect is destructive. It is noted that for fi+ there still are Cabibbo- 
suppressed WE diagrams, but for S+ + there are only PI diagrams. Therefore the predicted 
lifetime of S+ + is larger than that of other two baryons. We also employ other values for 
parameters a p , a\ and find that the resultant values can vary within 20% uncertainty. 

Our results are 

r(~+ c ) = 0.25 ps r(~++) = 0.67 ps and r(fi+) = 0.21 ps. 

These are generally consistent with the results obtained by Kiselev et al. j3] and Guberina 
et al.Jsj], even though they used different models for calculating the hadronic matrix elements. 
Concretely, they used the diquark picture and attributed the non-perturbative effects into 
the wavefunction of the diaquark at origin. Kiselev et al. gave t(H+) ~ 0.16 — 0.22 ps 
r(S++) ~ 0.40 - 0.65 ps and r(fi+) ~ 0.24 - 0.28. 

Although all the theoretical predictions based on different models agree with each other, 
they are obviously one order larger than the upper limit of the measured value on the 
lifetime of E+ (0.033 ps) by the SELEX co 
some authors, may come from experiments 



2n 



laboration[6|. This deviation, as suggested by 



So far the difference between theoretical 



predictions and experimental data may imply some unknown physics mechanisms which 
drastically change the value, if the future experiment, say at LHCb, confirms the measure- 
ment of the SELEX. Recently, several groups have studied the possibility of doubly heavy 



baryon production at hadron collider LHC and future linear collider ILC|29|, 301] and the 
effective field theories for two heavy quarks system are also further investigated 3l| . We are 
expecting the new data from more accurate experiments at LHC and ILC to improve our 
theoretical framework and determine if there are contributions from new physics beyond 
the standard model. 
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FIG. 1: non-spectator effects contribution to lifetime of ^ C cd 
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FIG. 2: non-spectator effects contribution to lifetime of 
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FIG. 3: non-spectator effects contribution to lifetime of Q, 
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